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Abstract 



\^ We study the polarization-polarization and polarization-temperature corre- 

lations in standard adiabatic scenarios for structure formation. Temperature 
anisotropies due to gravitational potential wells and oscillations in the photon- 

o 

baryon-electron fluid on the surface of last scattering are each associated with 
a correlated polarization pattern. While the 'correlated part' of the polariza- 
tion has an r.m.s. of only a third of the total signal, it may still be measurable 

I 

O by mapping a large area on the sky. We calculate the expected signal to noise 

^ ratio for various measures of the polarization in a hypothetical mapping ex- 

periment such as those now being planned. 



I. INTRODUCTION 

Since COBE's detection of anisotropies in the cosmic background radiation (CBR) [1], 
attention has focussed on obtaining higher resolution measurements, and ultimately maps 
of the temperature anisotropies on the sky. The angular correlations and statistics of the 
anisotropy pattern will yield valuable clues as to the mechanism or mechanisms involved 
in the formation of large scale structure in the universe - whether the perturbations were 
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adiabatic or isocurvature, gaussian or nongaussian, whether there was a significant gravity 
wave component, or cosmic defects were involved. However there is substantial degeneracy 
amongst the predictions of different theories [2], [3], [4], and it is worth asking whether any 
additional information that might further discriminate between theories could be extracted 
from the microwave sky. 

The idea that the polarization of the microwave sky might provide such additional in- 
formation is not new (see e.g. [5], [6], [7], [8]), and some experimental limits have already 
been set [9], [fO], [ff]. The expected level of linear polarization is low (typically 5% of the 
anisotropy), but with experiments currently being planned to map the sky temperature to 
an accuracy of 3/iK per pixel one is clearly close to the level required for a measurement. 
Despite the low amplitude, it has long been realized that there is a compensating factor for 
ground-based experiments, namely that in principle one has only to measure a single point 
on the sky, and in this case both polarizations are affected similarly by the atmosphere. 

In this Letter, we extend previous theoretical work [7] to the temperature-polarization 
correlation function in standard adiabatic = 1 scenarios. We discuss the advantages the 
temperature-polarization correlation has in situations where the detector noise per pixel 
is larger than the signal, and show that temperature-polarization cross-correlation may be 
measurable in experiments with large sky coverage such as are currently being proposed. 

Why should the polarization of the sky be correlated with the temperature anisotropy? 
At redshifts > 1300, the photons, baryons and electrons constitute a tightly coupled fiuid. 
Acoustic oscillations of this fiuid produced as perturbation modes cross the horizon lead to 
the 'intrinsic' and 'Doppler' temperature fiuctuations which dominate the CBR anisotropics 
on small angular scales. The linear polarization produced in radiation which Thomson scat- 
ters off an electron is proportional to the quadrupole moment of the incident photon phase 
space density. This is vanishingly small in the fiuid limit. But as photons decouple, re- 
gions in which there is a converging or diverging velocity field develop a quadrupole of order 
rSr ~ — rVur, with r the mean free time, 6^ and the density and velocity perturbation. 
A positive fiuctuation in the radiation density 6^ produces a corresponding fiuctuation in 
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the temperature on the sky (hotter regions decouple later, and photons from them redshift 
less). Depending upon whether 6^ is positive or negative, such a region may have a converg- 
ing or a diverging velocity field. Around a temperature hot spot one would find a radial 
polarization pattern if 6^ were positive, tangential if it were negative. Note that all modes 
of a given wavelength oscillate in phase, so different signs of the temperature-polarization 
cross correlation are expected on different angular scales. 
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FIG. 1. A gravitational potential hill on the surface of last scattering leads to a hot spot on 
the sky since photons (71, 72) blueshift as they leave it. Blueshifted photons scattered to us from 
one side, such as 72, would be linearly polarized into the page, creating a tangential pattern about 
the hot spot. 

On larger angular scales, CBR anisotropics due to gravitational potential fiuctuations 
on the surface of last scattering dominate in the simplest adiabatic scenarios. These also 
produce a correlated polarization signal (Figure f), a tangential pattern of polarization 
around temperature hot spots, radial around cold spots. The scale of the polarization 
pattern is limited by the photon mean free path at last scattering, of order the horizon at 
that time - which subtends an angle of 2° with standard recombination, and 6° for a fully 
ionized universe. 

II. CALCULATION 

To compute the two point correlation functions, we evolve the photon distribution func- 
tion f(x,p,t) [12]. We assume the initial perturbations are adiabatic, pure growing modes. 
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with the matter and radiation being described by perfect fluids. Here we shall present re- 
sults for an = 1, 0^ = 0.05, h = 0.5 CDM dominated universe with standard thermal 
history, and a universe with no recombination (an extreme example of a reionized universe). 
Generalizations to other cosmological parameters are straightforward and will be presented 
elsewhere [13]. 

The initially Planckian, unpolarized photon distribution function is evolved forward using 
the general relativistic Boltzmann equation for radiative transfer with a Thomson scattering 
source term [14]. For simplicity we consider only scalar perturbations. In this case one needs 
to evolve only two transfer equations, for the components of f corresponding to Stokes' 
parameters / and Q [15]. These equations are written in terms of brightness functions, [16] 
A, ^ iSf,/ (ToI^), where To is the mean CMB temperature, / is the unperturbed Planck 
distribution, and i = T^P . 

The computational scheme is that of Bond and Efstathiou [7], [12]. We expand the cos- 
mological perturbations in plane waves, and the brightness functions in spherical harmonics, 
converting the transfer equations to a hierarchy of ordinary differential equations. 

Once evolved to the present epoch, the Legendre expansion coefficients A) (z = T, P) are 
used to construct the temperature and polarization correlation functions 







l>2 
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where (^, cf)) are the usual spherical polar angles, q = (sin ^ cos ^, sin ^ sin ^, cos ^), the axes 



used to define the Stokes parameters are e^; and e^^, and 




(6) 
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cr = jedk\Ar\' (7) 

TP ~ 



{1 + 2] 



^pp _ - 4) 



(^ + 4) 



J2i2l' + 1) 1 edkA^Af*au' (8) 
^(2/' + l) / edkAf,Af*hi,. (9) 

The constants a;;/ and a;;/ are given by a;;/ = j\ dxPi[x)Pj^[x)^ and a;;/ = j\ dxPi{x)P^{x) 
which have simple closed form expressions. In deriving equations (4) and (5) we further 
assumed ^ ^ 1. 

Figure 2 shows the QT cross correlation function for ^ = 0. The cos 2cf) dependence of 
{QT) means that it must vanish at zero angular separation in order to be single valued. Scales 
< 2° are subtended by physical wavelengths shorter than the horizon at recombination, where 
the oscillation of the baryon and radiation fluids is reflected in the structure of the two 
point function. At larger angular scales, the negative tail expected from the 'gravitational 
potential' argument above is seen. 




FIG. 2. QT correlation function in units of (Txctq, where ctj is the temperature anisotropy 
variance and ctq is the variance of the total polarization. 

If the primordial perturbations are Gaussian, so are the temperature and polarization 
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fields. In this case tlie statistics are completely described by these two point functions. 
We have constructed realizations of the microwave sky by making use of the small angle 
approximation, in which spherical harmonics may be replaced by plane waves, 

r(6/) = ^r(n)e''2"^ • (10) 
n 

and similarly for Q and [/, where is the angular size of the map, and the wavenumber 
n has integer components. Each Fourier mode n of the temperature anisotropy field T 
is allocated a random Gaussian-distributed, complex variable (i with zero mean, and unit 
variance (CiCi) =1? according to T[n) = \J Cf(i , where / ~ (27r/n0). A second independent 
random variable 1^2 is used to define the Q Fourier modes: 

Q{n) = (gc(n)Ci + Qu{n)C2) cos 2(^„ (11) 

where Qc{n) and Qu{n) are chosen to reproduce the {QT) and {QQ) correlation functions. 
The U modes are allocated in a similar way. 

The polarization field is therefore a sum of two components, which are respectively 
correlated, Qc{0), and uncorrelated, Qjj{6), with the temperature anisotropy. Fig 3 shows 
the correlated component overlaid on the temperature field. The length of each vector is 
proportional to [Ql;{6)-\-Ul;{6)\^ and the orientation is given by 2^ = {U c I Q c) ■ There 
are clear correspondences between the temperature and the correlated component of the 
polarization. The most obvious features are seen around hot and cold spots in an otherwise 
uniform background, with radial and tangential polarization patterns respectively. Note that 
the Figure only shows the correlated component of the polarization - the effect of including 
the larger uncorrelated component is to mask almost entirely all obvious correspondences 
with the temperature map. 

III. INTERPRETATION 

The correlated part of the polarization is a small part of the total signal. Figure 4 
shows the power spectra of the temperature-correlated polarization, {QcQc)^ and the total 
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FIG. 3. 2° X 2° temperature map (smoothed with a Gaussian beam with FWHM of 10') with the 
correlated component of the polarization overlaid. Note the cold spot at (95,35) with a surrounding 
tangential polarization pattern, and the hot spot at (25,65) with a radial pattern. 



7 



polarization, {QQ). The variance of the total polarization, (jg, corresponding to the area 
beneath the curve, is approximately seven times that of its correlated component, Cg^. (The 
ratio for the constantly ionized universe is similar.) From a measured map of the temperature 
anisotropy, it is straightforward to construct a map of the correlated part of the polarization, 
as we have done in Figure 2. Because the uncorrelated part of the polarization is large, the 
map of the correlated part is a useful predictor of the total polarization only in a statistical 
sense. A one to one correspondence between features in the temperature-correlated map 
and those in the total polarization is not expected. However, one is much more likely to find 
a peak in the total polarization at a peak in the correlated map than at a random point in 
the sky. 




FIG. 4. The power spectra of the total and correlated polarization in a universe with a stan- 
dard thermal history normalized to the COBE measurement of the temperature variance with 10° 
smoothing. 

This statement can be quantified. The probability of finding an tktq peak in the total 
polarization if one is at an maq^ peak in the correlated polarization is 

P = , 2 ^ 2 ,,1 exp[-(n(7Q - maQ^f/2{al - a^J], (12) 

(27r((7^ -(7^J)2 
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while the probability of finding an na peak by looking at a random point in the sky is 
e~"^/^/(27r(Tg)2 . Thus, for example, the odds of finding a Ictq peak in the total polarization 
at a 3(TQp peak in the correlated polarization are three times greater than at a random point. 
If degree scale temperature anisotropy maps become available, it may be useful construct 
the map of correlated polarization (in a given theoretical scenario), and use it as a guide to 
the 'best' points at which to observe the polarization. 

Although one expects the correlated polarization to be small, one can attempt to measure 
the {QT) correlation function directly. In fact, measuring {QT) could prove easier than 
measuring {Q'^) because it is less susceptible to noise in the polarization measurement. 
Consider a detection obtained from N measurements of the polarization, Qi it (7/5, where 
cr£) is the detector noise. If the measurements are sufficiently isolated from each other to be 
uncorrelated, then the measured variance will be, 

Q'^lmeas = ± ^l^Wl + Crh] (13) 

where ctq is the true polarization variance. Whereas, the measured temperature-polarization 
correlation will be 

/ 2 

QT{9,^)U,, = {QT{9,^)) ±^^[al + al]^ (14) 

where cr^ is the variance of the temperature anisotropy. (We have assumed that the noise 
in the temperature anisotropy is negligible.) In the limit of large detector noise, the error 
in measuring grows as ajj while the error in QT grows as (7/5, so that it becomes more 
advantageous to search for QT correlations. 

If one has a full sky temperature map but only a sparsely sampled polarization map, 
the noise in the temperature-polarization correlation can be reduced by including all of the 
temperature information. For example, we can calculate the expected correlation between 
the polarization at a given point and the temperature in a ring of radius 6 about that point. 
Since the QT correlation is proportional to cos 2^, if we define T{9) = ^ / d(/)T{9, (/)) cos 2^, 
then {QT[6)) = {QT[6^ cf) = 0°))/2. The error in measuring {QT[6)) is reduced significantly 
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because T[6) is an average over a number of uncorrelated patches. Thus, the variance in 
T{9), {T^{9)) = f d<f>cos2<f>CT{9{2 — 2cos^)), is much smaller than the variance in the 
temperature anisotropy itself. For a ring of radius ^, T is a weighted average of / ~ irO/Oc 
independent patches (where 6c is the temperature correlation angle) and its variance is 
(T'^{6)) ~ (Ty/2/. Assuming the rings about the polarization measurements do not have 
substantial overlap, this increases the signal to noise of {QT[9)) by a factor of ^//2 over 
that of {QT{e)). 

Alternatively, one can test the correlated polarization map (constructed as described 
above) by measuring the weighted average, {Q^q^°^^Q). For N uncorrelated measurements 
of the polarization. 

Comparing the signal to noise of {Q"^) and {Q^c^°^^Q)^ and using ~ ctq/T, one finds that 
if the noise is greater than the polarization signal, cro > 1.5(Tq, then it becomes easier to 
measure {QcQ) than {Q'^). 

Consider a hypothetical experiment which measures the polarization with a 0.5° FWHM 
beam of 1000 well separated patches on the sky. The expected polarization for a standard 
recombination model is ctq ~ 1.4/iA^, while the correlated polarization signal is about a 
third of this. If the noise level can be reduced to 3 fiK per pixel, the signal to noise of the 
polarization variance, (g2),and that of {Q'^'"''" Q) are both about five to one, while the signal 
to noise of the {QT[6)) is approximately three to one on an angular scale 6 ~ 1.3°. For a 
fully sampled polarization map the signal to noise would improve. For the no recombination 
case, the expected noise levels are comparable to these because even though the variance 
of the polarization on 0.5° is larger (~ 3/iA^), the correlation angle is also larger and fewer 
independent measurements can be made from the same sky coverage. These numbers ignore 
any other non-cosmological sources of polarization, but on such large angular scales one can 
hope that astrophysical sources e.g. radio galaxies, hot gas in clusters, do not contribute 
significantly. 
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Given a fixed integration time, the optimal observation strategies for measuring the 
temperature polarization correlation and the polarization auto-correlation differ significantly. 
If the noise cro is statistical, i.e. inversely proportional to the square root of the time spent 
in making the observation, then the best strategy for measuring the polarization variance is 
to spend enough time at each observation such that cro is comparable to ctq before moving 
on. The optimal strategy for measuring {QT[6))^ however, is to maximize the number of 
observations, independent of the noise level, i.e. to map out the polarization over the whole 
sky. Such a strategy might most easily be realized in conjunction with a full sky anisotropy 
mapping experiment such as those presently being planned. 

To conclude, the temperature-polarization correlation provides another observable quan- 
tity which may be used to probe the physics of the density perturbations on the surface 
of last scattering. Although the magnitude of the effect is small, we have shown that a 
reasonable signal to noise ratio is within reach of projected experiments. It would be very 
interesting to extend the calculations reported here to other structure formation scenarios, 
such as baryon isocurvature and defect models. 

We thank K. Ganga, P.J.E. Peebles and D. Wilkinson for useful conversations. RC 
thanks R. Davis and P. Steinhardt for their collaboration in the initial development of the 
photon evolution code. The work of DC was supported by EPSRC (UK), while that of RC 
and NT was partially supported by NSF contract PIIY90-21984, and the David and Lucile 
Packard Foundation. 
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